Relevance of electron-lattice coupling in cuprates superconductors 
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We study the oxygen isotope ( 16 0, 18 0) and finite size effects in Yi- X Pr^I^C^Or-a by in- 
plane penetration depth (X a b) measurements. A significant change of the length L c of the su- 
perconducting domains along the c-axis and \\ b is deduced, yielding the relative isotope shift 
AL C /Lc ~ AXa b /Xa b ~ —0.14 for x — 0, 0.2 and 0.3. This uncovers the existence and relevance of 
the coupling between the superfluid, lattice distortions and anharmonic phonons which involve the 
oxygen lattice degrees of freedom. 

PACS numbers: 74.72.Bk, 74.25.Kc, 74.62.Yb 
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Since the discovery of superconductivity in cuprates by 
Bcdnorz and Muller|lJ a tremendous amount of work has 
been devoted to their characterization. The issue of inho- 
mogeneity and its characterization is essential for several 
reasons, including: First, if inhomogeneity is an intrin- 
sic property, a re-interpretation of experiments, measur- 
ing an average of the electronic properties, is unavoid- 
able. Second, inhomogeneity may point to a microscopic 
phase separation, i.e. superconducting domains, em- 
bedded in a nonsuperconducting matrix. Third, there 
is neutron spectroscopic evidence for nanoscale cluster 
formation and percolative superconductivity in various 
cuprates@>Q. Fourth, nanoscale spatial variations in the 
electronic characteristics have been observed in under- 
doped Bi2Sr2CaCu90»+*_ with scanning tunnelling mi- 
croscopy (STM)0, la, la 0- They reveal a spatial segre- 
gation of the electronic structure into 3nm diameter su- 
perconducting domains in an electronically distinct back- 
ground. On the contrary, a large degree of homogeneity 
has been observed by Renner and Fischer [8J. As STM is a 
surface probe the relevance of these observations for bulk 
and thermodynamic properties remains to be clarified. 
Fifth, in YBa 2 Cu 3 07-5, MgB 2 , 2H-NbSe 2 and Nb 77 Zr 23 
considerably larger domains have been established. The 
magnetic field induced finite size effect revealed lower 
bounds ranging from L = 182 to 814A (3, OH- Sixth, 
since the change of the lattice parameters upon oxygen 
isotope exchange is negligibly small [ill fl2j| . the occur- 
rence of a significant modification of the domains spatial 
extent, will provide clear evidence for superconductivity 
mediated by local lattice distortions. 

This letter concentrates on Yi- X Pr.,;Ba 2 Cu307_,5 with 
16 and 18 0. It addresses these issues by providing clear 
evidence for a finite size effect on the in-plane London 
penetration depth, revealing the existence of supercon- 
ducting domains with spatial nanoscale extent and its 
significant change upon oxygen isotope exchange. 

Polycrystalline samples of Yi_ 3; Pr 2 ;Ba 2 Cu307_5 (x = 
0, 0.2, 0.3) were synthesized by solid-state reactions 
and their phase-purity was examined using powder x- 
ray diffraction. For each doping concentration the sam- 



ples were reground in a mortar for about 60min. Pow- 
der samples with a grain size of < 10/im were obtained 
by using a system of 20/15/10^m sieves. Oxygen iso- 
tope exchange was performed by heating the powder 
in 18 2 gas. In order to ensure the same thermal his- 
tory of the substituted ( 18 0) and not substituted ( 16 0) 
samples, two experiments (in 16 Q 2 and 18 2 ) were al- 
ways performed simultaneously |13j. To achieve com- 
plete oxidation the exchange processes were carried out 
at 550°C during 30 h, followed by slow cooling (20°C/h). 
The 18 O content, determined from a change of the sam- 
ple weight after the isotope exchange, was found to be 
89(2)% for all samples. Field-cooled (FC) magnetization 
measurements were performed with a Quantum Design 
SQUID magnetometer in field range 0.5 to lOmT and 
a temperature range 5K to 100K. The powder samples 
(~100mg) were put in a quartz ampule. To guarantee 
the same experimental conditions (sample geometry and 
the background signal from the sample holder) , the same 
ampule was used. The absence of weak links between 
grains was confirmed by the linear scaling of the FC mag- 
netization measured at 5K in 0.5mT, lmT and 1.5mT. 
The Meissner fraction / was calculated from the mass 
and the x-ray density, assuming spherical grains. An ex- 
ample of / versus T is displayed in the insert of Fig^ 
Assuming spherical grains of radius R, the data were an- 
alyzed on the basis of the Shoenberg formula 0| , allow- 
ing to calculate the temperature dependence of the effec- 
tive penetration depth X e ff (T) /Xeff (0). For sufficiently 
anisotropic extreme type II superconductors, including 
Yi_ 2; Pr 2 ;Ba 2 Cu307_ ( 5 , X e ff is proportional to the in- 
plane penetration depth, so that A e // = 1.31A„?,[lfil llfij. 
Since it is not possible to extract the absolute value of 
X a b from our measurements we take X a b (Oif) from /zSR 
measurements 1 1 7| to normalize the data. In the main 
panel of Fig^ we displayed the resulting temperature 
dependence of 16 A 2 b (0) /X 2 ab (T) for the 16 and 18 
samples of Yo yPro 3Ba 2 Cu307_a. For comparison we 
included 16 Xl b (0) /A 2 b (T) = 16 A 2 b (0) /A 2 , afc (1 - T/T c f 
with the critical exponent v = 2/3, 16 T C = 5A.9K and 
18 T C = 53.7K to indicate the asymptotic critical behav- 
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FIG. 1: ( 16 A a6 (0)/A a6 (T)) 2 versus T for the 16 and 18 
samples of Yo.7Pro.3Ba2Cu307_a. The solid lines indicate 
the leading critical behavior of the homogeneous systems as 
explained in the text. Insert: Meissner fraction / versus tem- 
perature (FC lmT). The error bars are smaller than the size 
of the data points. 



ior for an infinite superconducting domain. Apparently, 
the data are inconsistent with such a sharp transition. It 
clearly uncovers a rounded phase transition which occurs 
smoothly and with that a finite size effect at work. In this 
context it is important to emphasize that this finite size 
effect is not an artefact of Yi_a;Pr x Ba2Cu307_,5 or the 
particular technique used to evaluate 1/A^ b (T). Indeed, 
this rounding is also seen in the data for YBa2Cu30y 
llit. La 2 _ a; Sr 2; Cu04 with x=0.1, 0.15 and O.20, and 
Bi 2 Sr2CaCu208+(5 single crystalsjitj HJ. Moreover, in- 
dependent evidence for superconducting domains of finite 
extent, stems from the analysis of specific heat data|l0|. 
Given the mounting evidence for these domains, their be- 
havior upon oxygen isotope exchange is expected to offer 
valuable clues on the relevance of local lattice distortions 
in the mechanism mediating superconductivity. 

To elucidate this issue we perform a finite size scal- 
ing analysis of the in-plane penetration depth data 
for Yi_ 2; Pr 2 ;Ba2Cu307_ ( 5 with 16 and 18 0. Supposing 
that cuprate superconductors are granular, consisting of 
spatial superconducting domains, embedded in a non- 
superconducting matrix and with spatial extent L a , L b 
and L c along the crystallographic a, b and c-axes, the cor- 
relation length £j in direction i, increasing strongly when 
T c is approached cannot grow beyond Li. Consequently, 
for finite superconducting domains, the thermodynamic 
quantities like the specific heat and penetration depth 
are smooth functions of temperature. As a remnant of 
the singularity at T c these quantities exhibit a so called 
finite size effect 22], namely a maximum or an inflection 



point at T Pi , where ^ (T Pi ) — Li. There is mounting ex- 
perimental evidence that for the accessible temperature 
ranges, the effective finite temperature critical behavior 
of the cuprates is controlled by the critical point of un- 
charged superfluids f3D-XY)[iJ. l21j. In this case there is 
the universal relationship 
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between the London penetration depth Aj and the trans- 
verse correlation length £| in direction i(3,H^|. As afore- 
mentioned, when the superconductor is inhomogeneous, 
consisting of superconducting domains with length scales 
Li, embedded in a non-superconducting matrix, the £*'s 
do not diverge but are bounded by 
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A characteristic feature of the resulting finite size effect 
is the occurrence of an inflection point at T Pk below T c , 
the transition temperature of the homogeneous system. 
Here 



and Eq.0J reduces to 



A, (T)Xj(T) 
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In the homogeneous case 1/ (Aj (T) Xj (T)) decreases con- 
tinuously with increasing temperature and vanishes at T c , 
while for superconducting domains, embedded in a non- 
superconducting matrix, it does no vanish and exhibits 
an inflection point at T Pk < T c , so that 
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We are now prepared to perform the finite 
size scaling analysis of the penetration depth 
data. In FigH we displayed w X 2 ab (T = 0) /X 2 ab (T) 
and d(Xl b (T = 0)/Xl b (T))/dT versus T for 
Yo.7Pro.3Ba 2 Cu 3 07_5. The solid lines are 
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v 3 A Qb (0)/ lb A Qb( 

( 16 A ab (0)/ 18 A Qb (T)) z 1.4(1 
v = 2/3, 16 T C = 54.9K, 18 T C = 53. 7K and the 
dash dot lines the corresponding derivatives, indicating 
the leading critical behavior of a domain, infinite in 
the c-direction. The extreme in the first derivative 
around T « 52. IK and 51K for 16 and 18 respec- 
tively clearly reveal the existence of an inflection point, 
characterizing the occurrence of a finite size effect in 
l/A^ 6 (Eq.0). Using Eq.(@J and the estimates for 
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A^(T P J/ 16 A^(0), 18 A2 b (T p J/^ b(0 ) and 
A ab (0) listed in Table I, we obtain 16 L c = 19.5 (8)A 
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FIG. 2: 16 \ 2 ab (T = 0)/\ 2 ab (T) (Q : 16 0, . : 18 0) and 
d (X 2 ab (T = 0) /\ 2 ab (T)) /dT (□ : 16 O, ■ : 18 0)versus T for 
Yo. 7 Pro.3Ba 2 Cu3 7 with 16 and 18 0. The solid and dash 
dot lines indicate the leading critical behavior of the homoge- 
neous system as explained in the text. 
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FIG. 3: 16 X 2 ab (T = 0)/X 2 ab (T) (O 
d (X 2 ab (T = 0) /\ 2 ab (T)) /dT (□ : 16 O, ■ : 18 0) versus T for 
YBa 2 Cu 3 07 with ie O and ls O. The solid and dash dot lines 
indicate the leading critical behavior of the homogeneous sys- 
tem as explained in the text. The dot lines are quadratic fits 
used to determine the inflection points T Pc . 



and x — 0.2 and extracted the in-plane penetra- 
tion depth, as outlined above. In Fig|3| we show 
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K b {T = 0)/Xl b {T) and d (X 2 ab (T = Q) / X 2 b (T)) / dT 



versus T for YBa 2 Cu 3 7 _ 5 with 16 and 18 0. The solid 

' T/ 16 T C ) 2/3 , 

,2/3 



lines are 



( 16 A ab (0)/ 16 A Qb (T))^ = 3.9(1 

= 3.55 (l - T/ 18 T C 



( 16 A ab (0)/ 18 A ab (T)) z = 3.55(1-T/^T C )"^ with 
W T C « 18 T c = 90.3K and the dash dot lines the cor- 
responding derivatives, indicating the leading critical be- 
havior of the homogeneous system. The finite size esti- 
mates for various quantities and the corresponding iso- 
tope shifts, for an 18 content 89% are summarized in 
Table I (we define the relative oxygen isotope shift of 
a physical quantity X as AX/X = ( 18 A - 16 X)/ W X). 
While the isotope effect on T c and the inflection points 
T Pc is very small, there is an appreciable shift of l/X 2 b . 
As indicated in FigEI we used a quadratic fit around the 
inflection point to determine T Pc . Nevertheless, the error 
left (0.1 to 0.2K) leads to the main uncertainty in terms 
of A 2 b (T Pc ) /X 2 b (0). From Table I several observations 
emerge. First, L c increases systematically with reduced 
T Pc . Second, L c grows with increasing x and upon iso- 
tope exchange ( 16 0, 18 0). Third, the relative shift of T Pa 
is very small. This reflects the fact that the change of 
L c is essentially due to the superfluid, probed in terms 

AA 2 b /A 2 b for x = 0, 0.2 
2 1 and 



of A 2 b . Accordingly, AL C /L. 



Indeed the relative shifts of T Pc , A 2 b (T Pc 



and 0.3 

L c are not independent. Eq.Q implies 
AL C _ AT Pc AX 2 ab (T Pe 



K b (T P J 
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X 





0.2 


0.3 


AT Pc /T Pc 


-0.000(2) 


-0.015(3) 


-0.021(5) 




0.12(5) 


0.13(6) 


0.16(5) 


AA^ b (T p J/A- b (T p J 


0.11(5) 


0.15(6) 


0.15(5) 




4.4(2) 


4.0(2) 


4.4(2) 




4.9(2) 


4.6(2) 


5.2(2) 


W T Pc (K) 


89.0(1) 


67.0(1) 


52.1(2) 


l "T Pc (K) 


89.0(1) 


66.0(2) 


51.0(2) 


(A) 


9.7(4) 


14.2(7) 


19.5(8) 


(A) 


10.9(4) 


16.0(7) 


22.6(9) 


16 A Qb (0)(A) 


1250(10) 


1820(20) 


2310(30) 



and 18 L C = 22.6 (9)A for the spatial extent of the 
superconducting domains along the c-axis. Note that 
the rather broad peak around the inflection point reflects 
the small value of L c . Indeed, in Bi2212, where the 
same analysis gives L r ~ 68A, this peak was found to be 
considerably sharper [24|. 

To explore the dependence of this change on Pr 
concentration we performed analogous magnetization 
measurements on Yi_ 2; Pr 2 ;Ba2Cu307_ ( 5 with x = 



Table I: Finite size estimates for 16 T Pc , 16 T Pc , 
AT Pc T Pc - 18 T Pc , 16 A 2 b ( 16 T Pc )/ 16 A 2 b (0) and 

18 Aab ( 18t pc) / 16 A^ b (0), and the resulting relative shifts 
AT Pc / 16 T Pc and AA 2 ,/ 1 ^ ( 16 T Pc ) for an 18 content 
of 89%. 16 L Pc , 18 L Pc and AL p J 16 L^ c follow then via 
Eq.(@}. 16 A ab (0) are /iSR estimates" 

To appreciate the implications of these estimates, we 
note that for fixed Pr concentration the lattice param- 
eters remain essentially unaffected [Til ll^. Accordingly, 
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an electronic mechanism, without coupling to local lattice 
distortions and anharmonic phonons, implies AL C = 0. 
On the contrary, a significant change of L Pc upon oxy- 
gen exchange uncovers the coupling to local lattice dis- 
tortions and anharmonic phonons involving the oxygen 
lattice degrees of freedom. A glance to Table I shows 
that the relative change of the domains along the c- 
axis upon oxygen isotope exchange is significant, ranging 
from 12 to 16%, while the relative change of the inflec- 
tion point or the transition temperature is an order of 
magnitude smaller. For this reason the significant rela- 
tive change of L c at fixed Pr concentration is accompa- 
nied by essentially the same relative change of \ 2 ab , which 
probes the superfluid. This uncovers unambiguously the 
existence and relevance of the coupling between the su- 
perfluid, lattice distortions and anharmonic phonons in- 
volving the oxygen lattice degrees of freedom. Potential 
candidates are the Cu-0 bond-stretching-type phonons 
showing temperature dependence, which parallels that 
of the superconductive order parameter}^. Independent 
evidence for the shrinkage of limiting length scales upon 
isotope exchange stems from the behavior close to the 
quantum superconductor to insulator transition where T c 
vanishes j!^]. Here the cuprates become essentially two 
dimensional and correspond to a stack of independent 
slabs of thickness d s [j?! Eij ] . It was found that the rela- 
tive shift Ad s /d s upon isotope exchan ge a dopts a rather 
unique value, namely Ad s /d s « —0.03 26] . 

In conclusion, we reported the first observation of 
the combined finite size and oxygen isotope exchange 
effects on the spatial extent L c of the superconduct- 
ing domains along the c-axis and the in-plane penetra- 
tion depth A a h. Although the majority opinion on the 
mechanism of superconductivity in the cuprates is that 
it occurs via a purely electronic mechanism involving 
spin excitations, and lattice degrees of freedom are sup- 
posed to be irrelevant, we have shown the relative isotope 
shift AL C / 16 L C w AX 2 ab / w X 2 ab w -0.15 uncovers clearly 
the existence and relevance of the coupling between the 
superfluid, lattice distortions and anharmonic phonons 
which involve the oxygen lattice degrees of freedom. 
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National Science Foundation. 
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